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Abstract 

We prove the existence of a countable family of globally regular solutions of spherically 
symmetric Einstein-Klein-Gordon equations. These solutions, known as mini-boson stars, 
were discovered numerically many years ago. 
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r <P : 1 Introduction 

Boson stars are compact gravitationally bound soliton-like equilibrium configurations of bosonic 
fields. The simplest kind of boson star, which is made up of a self-gravitating free complex 
massive scalar field, was conceived over thirty years ago by Kaup jy and Ruffini and Bonazzola [^) 
who found numerically the ground state solution to the spherically symmetric Einstein-Klein- 
' Gordon (EKG) equations. A decade later the systematic numerical analysis of these equations 

was performed by Friedberg, Lee, and Pang || who rediscovered and extended the results of jl[ [| , 
' in particular they found a countable sequence of excited states. 

The aim of this paper is to give a rigorous proof of existence of solutions found in ^J, |j] . In 
the physics literature these solutions are usually referred to as mini-boson stars ('mini' because 
they are tiny objects with mass ~ where m is the boson mass). What are the motivations 
for studying such objects? Let us mention three possible reasons varying from physical to purely 
I mathematical. First, most theories of elementary particles predict the existence of massive 

^Jy bosons which interact weakly with baryonic matter. To the extent one believes in these models, 

one should accept their consequences, like boson stars. From this standpoint, the recent surge 
of interest in boson stars is largely due to the suggestion that the dark matter could be bosonic 
since then some fraction of the missing mass of the universe would float around in the form 
of boson stars. Second, even if massive scalar fields do not exist in nature, they provide one 
of the simplest fundamental matter sources for the Einstein equations and, as such, are ideal 
theoretical "laboratories" for studying the dynamics of gravitational collapse. Mathematically, 
these studies amount to the analysis of the Cauchy problem for the EKG equations. Boson 
stars play an important role in this context as candidates for intermediate or final attractors 
of dynamical evolution. Finally, and admittedly most interestingly for us, mini-boson stars are 
non-perturbative solutions of the EKG equations in the sense that they have no regular flat- 
spacetime limit (one manifestation of this property is the fact mentioned above that the total 
mass of a mini-boson star is inversely proportional to the gravitational constant G). In this 
respect mini-boson stars are similar to the Bartnik-McKinnon solutions of the Einstein- Yang- 
Mills equations [[|. However, in contrast to the Bartnik-McKinnon solutions, the mini-boson 
stars are not static: although the metric and the stress-energy tensor of the scalar field are 
time-independent, the scalar field iself has the form of a standing wave <p(r,t) — e tuJt 4>(r). This 
fact has an important consequence at the ode level, namely the lapse function does not decouple 
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from the Klein-Gordon equation and the hamiltonian constraint which means that we have to 
deal with a 4-dimensional (nonautonomous) dynamical systemQ Below we analyze this system 
using a shooting method which is similar in spirit (but quite different in implementation) to the 
proof of existence of the Bartnik-McKinnon solutions || . 

The paper is organized as follows. In Section 2 we derive the field equations together with 
the boundary conditions and discuss some basic properties of solutions. We also formulate the 
main theorem and sketch the heuristic idea of its proof. In Section 3 we prove the local existence 
of solutions near the origin. In Sections 4 and 5 we discuss the limiting behavior of solutions 
for small and large values of the shooting parameter, respectively. In Section 6 we derive the 
asymptotics of globally regular solutions. Section 7 contains some technical results concerning 
the behavior of singular solutions. Finally, in Section 8, using the results of Sections 4-7, we 
complete the proof of the main theorem by a shooting argument. 



2 Preliminaries 

The action for the EKG system is given by 

S = J d^x^—g (-^ - ^d a pd a <t> - \m 2 r<^j , (2.1) 

where R is the scalar curvature of the spacetime metric g a b, 4> is the complex scalar field, and 
to is a real constant called the boson mass. The associated field equations are the Einstein 
equations 

Rab - \g ab R = 8nGT ab (2.2) 
with the stress-energy tensor of the scalar field 

Tab = \{d a ^d b ^ + datdbf) - ^g ab (g cd d c <f>*d d cl> + m 2 <^), (2.3) 
and the Klein-Gordon equation 

(□ - m 2 )(t) = 0, (2.4) 

where □ is the d'Alembertian operator associated with the metric g a b- Now, we assume that 
the fields are spherically symmetric. We write the metric using areal radial coordinate and polar 
slicing 

ds 2 = ~e- 25 Adt 2 + A- X dr 2 + r 2 dtl 2 , (2.5) 

where dQ, 2 is the standard metric on the unit two-sphere, and A and S are functions of (t,r). 
In this parametrization the (relevent components of) Einstein equations have the particularly 
simple form 

1 — A 

d r A = 87rGrT 00 , (2.6) 

r 

d r S = -iirGrA-^Too + Tn), (2.7) 

d t A = S-nGre~ s ATqi , (2.8) 

^^For comparison, the static spherically symmetric Einstein- Yang-Mills equations reduce (within the purely 
magnetic ansatz) to a 3-dimensional (nonautonomous) dynamical system. 
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where the components of stress-energy tensor T a b are expressed in the orthonormal frame deter- 
mined by the metric (2.5) (e = e s A~ 1 / 2 d t , e 1 = A 1 ' 2 d r ). From (2.3) we obtain 

Too = ^(Aldr^ + A-'e^ldt^+m 2 ^ 2 ), (2.9) 

Th = ^dr^ + A^e^dt^-mM 2 ), (2.10) 

Tqi - ^e s (d r <j>*d t <j> + d r <j>d t <j>*). (2.11) 

The remaining components of Einstein's equations are equivalent to the Klein-Gordon equation. 

For the scalar field cj) we assume the standing wave ansatz (f>(r, t) = cxp(iojt)<j)(r), where u is 
a real constant. Then, due to the U(l) symmetry of the action, the stress-energy tensor and the 
metric are time-independent. Morever, Til = so Eq. (2.8) is trivially satisfied. In terms of the 
dimcnsionless variables 

x = mr, f(x) = \Z47rG4>{r) (2.12) 

and the auxiliary variable 

C(x) = —A~ 1 e s , (2.13) 
m 

Eqs. (2. 4), (2. 6) and (2.7) reduce to the following system of ordinary differential equations (here- 
after prime denotes ^) 

^Pj = -AC 2 )/, (2T4a) 

A' = ^—A- x (Af' 2 + AC 2 f + f 2 ), (2.14b) 
x 

C = ^(A-l + x 2 ! 2 ). (2.14c) 

Instead of A, it is sometimes convenient to use the "mass" function M(x) defined by A(x) = 
1 - 2M(x)/x. From (2.14b) we have 

M{x) = ±J* s 2 {Af 2 + AC 2 ! 2 + f)ds. (2.15) 

A spacetime is said to be asymptotically flat if i5(oo) is finite and 

lim M(x) = M x < oo. (2.16) 

The limiting value is interpreted as the total mass of a solution (in our case it is measured 
in units <^;)- In Section 6 we will show that the finiteness of mass implies that C has a finite 
limit and / decays exponentially as x — > oo. 

Besides the singularity at infinity, the field equations (2.14) have the fixed singular point at 
x = and a moving singularity at x, where A(x) = 0. Regularity of solutions at x = requires 
the following behavior 

f(x) = a + 0(x 2 ), A{x) = l + 0(x 2 ), C(x) =a + 0(x 2 ), (2.17) 

where a = /(0) and a = C(0) are arbitrary parameters (assumed positive without loss of 
generality). In Sect. 3 we will show that these parameters determine uniquely a smooth local 
solution to Eqs. (2. 14). 
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Definition 2.1. The solution of Eqs.(2.14) starting at x = with the behavior (2.17) is called 
the a— orbit. 

In the following whenever we write 'a solution' we always mean the a-orbit. Also when we 
write that some property holds for all x we always mean for all x > 0. We will frequently refer 
to the behavior of a-orbits in the (/, /')-plane; when we write, say, that the a-orbit enters the 
first quadrant (Ql for brevity), we mean that the projection of the a-orbit in the (/, /')-plane 
does so. 

Definition 2.2. The a-orbits which exist for all x and are asymptotically flat are called globally 
regular. 

Now, we are ready to formulate our main result: 

Theorem 2.1. For each a > 1, there is a decreasing sequence of parameters a n (n — 0, 1, 2, ...) 
such that the corresponding a n -orbits are globally regular. The index n labels the number of nodes 
of the function f(x). 

This theorem makes rigorous the numerical results obtained in [|J ||, [| . Notice that although 
the a-orbits are determined by two parameters, only the parameter a has to be fine-tuned so the 
shooting is essentially one-dimensional. 

In order to prepare the ground for the proof of Theorem 2.1 we discuss now some elementary 
global properties of a-orbits. 

Lemma 2.2. A(x) < 1 for all x > unless f{x) = 0. 

Proof. From (2.14b), A'(xq) < if A(xq) = 1 so A cannot cross 1 from below. Since A(x) < 1 
for small x, the lemma follows. □ 



Lemma 2.3. An a-orbit exists as long as A(x) > 0. 

Proof. If A(x) > for x < x < oo, then \im x ^ s M(x) exists (because M' > and M(x) < 
x < x) so lini^^g A exists as well. We will show that the orbit can be continued beyond x 
provided that \im x ^x A(x) > 0. Since < A < 1, the only obstruction to extending the solution 
is the possibility that C,/, or /' might be unbounded. To see that / is bounded we note 
that (xA)' < 1 - x 2 Af' 2 . Choose e > such that A(x) > A(x)/2 for x — e < x < x; then 
(xA)' < 1 — (x — e) 2 A(x)f' 2 /2 and integrating from x — e < x gives that J^_ e f'(x) 2 < oo and 
hence, by the Cauchy-Schwarz inequality, J-_ \f'(x)\ < oo. Thus, / is bounded. This implies 
by Eq.(2.14c) that (InC)' is also bounded so both C and 1/C are bounded. Now, (2.14a) says 
that x 2 f'/C is bounded so /' is bounded. □ 

Remark. It follows from Lemma 2.3 that the only possible obstruction to extendability of a- 
orbits to arbitrarily large x is lim x ^ 2 A = for some x. If that happens we will say that the 
solution crashes at x. 

Let us define the function g — 1 — AC 2 . The following two properties of this function will play 
an important role in our discussion. 

Lemma 2.4. We have 

(a) Ifg(x) < 0, then g'(x) > 0; 

(b) If g(xo) > 0, then g(x) > for all x > xq. 
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Proof. A simple calculation yields 

g' = C 2 (±^+xAf' 2 -xgfy (2.18) 

The part (a) follows immediately from (2.18). To prove the part (b) note that g'(xi) > if 
g{x\) — 0, so g cannot cross zero from above. □ 

The restriction a > 1 in Theorem 2.1 can be easily seen as follows. Suppose that there is a 
globally regular solution with a < 1. Since g(0) = 1 — a 2 , it follows from Lemma 2.4 that g(x) is 
positive for all x. Multiplying Eq.(2.14a) by / and integrating by parts we get that //' > for 
all x, hence f 2 is monotone increasing which is obviously impossible for globally regular solutions 
(in fact such solutions crash at finite x as follows easily from Eq.(2.14b)). Thus we have 

Lemma 2.5. There are no (nontrivial) globally regular solutions for a < 1. 

Note that Lemma 2.5 implies in particular that there are no static (a — 0) globally regular 
solutions. In view of Lemma 2.5 from now on we always assume that a > 1. 

Definition 2.3. The rotation function 9(x, a) of an a-orbit is defined by 8(0, a) = 0, tan 9(x, a) = 
—f'(x)/f(x) and 9(x, a) is continuous in x. We will drop the second argument of 8 if there is no 
danger of confusion. 

Now we list the basic properties of the rotation function of a-orbits which we will need below. 
Lemma 2.6. For any nonnegative integer n we have: 

(a) If 8(x\) > (n + l/2)n for some x\, then 8(x) > (n + l/2)n for all x > x\. 

(b) If 8(x\) < nn for some Xi and g(x\) > 0, then 8(x) < nir for all x > x\. 

(c) There are at most two values of x with 8(x) = nn. 

Proof, (a) We note that 8'(x) = (f 2 - ff")/{f 2 + f' 2 ), so 8'{x) = 1 if 8(x) = (n + 1/2)tt. 

(b) If x > x\ and g(x\) > 0, then g(x) > by Lemma 2.4. Next, we note that 8'(x) = 
-g(x)/A(x) < if 8(x) = nn and g{x) > 0. If g{x) = then 6'{x) = but 8"{x) = 
-g'(x)/A(x) < since g'(x) > when g(x) = by Eq.(2.18). 

(c) The function 9(x) — nn changes sign at each zero for which g{x) ^ 0. From Lemma 2.4, g 
changes sign at most once. Thus, for n > 0, 8(0) — nn<0 and at x\, the first zero of 8(x) — nn, 
if g{xi) > then by part b) 8(x) — nn < for all x > x\. If g{xi) < then 8(x) — nn changes 
sign at Xi, and hence, at X2, the next zero of 8(x) — nn, g{x2) > and hence 8{x) — nn < for 
all x > x 2 . For n — 0, 8(0) — nn = 0, 8(x) > near x = and if 8(x\) = then g(xi) > 0, 
hence, 8(x) < for all x > x\. □ 

Before going into details of Sections 3-8, let us outline the main idea of the proof of Theo- 
rem 2.1. According to this theorem there exists a countable family of globally regular solutions 
distinguished by nodal class. We first show (section 3) that there is a continuous one-parameter 
family of local solutions depending on a = f(0); we all these solutions a-orbits. In Section 6 we 
show that an a-orbit that has bounded rotation and that is defined for all £ is a globally regular 
solution, that is, it has the correct asymptotic behavior as x — > oo. The existence of a-orbits 
with bounded rotation that are defined for all x is proven in each nodal class by an inductive 
application of a shooting argument. The zeroeth solution we construct has 8(x, ao) < n/2 for all 
x; the first solution has 8(x,a\) < in/2 (and greater than n/2 for large x), etc. This is shown 
in Fig. 1 . The crucial step of our argument is the control of behavior of a-orbits for large and 
small values of the parameter a. In Section 4 we show that for sufficiently small a the a-orbit 
has arbitrarily large rotation; more precisely, there is a number b n such that 8{x, a) > nn for 
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some x if a < b n . In contrast, we show in Section 5 that for a >> 1 the a-orbit exits Q4 directly 
to Ql (see Fig. 1). 

Now, to prove the existence of a globally regular solution in the zeroeth nodal class we let 
ao = inf{a| 9(x, a) < n/2 for all x for which the a-orbit is defined}. Note that ao > b\ > 0. We 
then prove that the ao-orbit is the globally regular solutions in the zeroeth nodal class. It is clear 
that the ao-orbit has rotation 6{x, ao) < 7r/2 for otherwise all nearby orbits would have rotation 
> it/2 which contradicts the definition of ao- It is also easy to see that the ao-orbit cannot exit 
Q4 to Ql because again, nearby orbits would also do so which contradicts the definition of a . 
Hence, the ao-orbit must stay in Q4; it either crashes or is defined for all x and is a globally 
regular solutions in the zeroeth nodal class. Thus, it remains to show that the ao-orbit does not 
crash. The (technical) crash lemma of Section 7 shows that if an orbit crashes in Q4 then nearby 
orbits either crash in Q4 or exit Q4 to Ql. Thus the ao-orbit cannot crash because nearby orbits 
would all be in {a\9(x, a) < n/2 for all x for which the a-orbit is defined} and ao would not be 
the infimum of that set. 

To show the existence of globally regular solutions in higher nodal classes we proceed as 
above. We let a n = mf{a\9(x,a) < (n + l/2)n for all x for which the a-orbit is defined}. We 
then show that 9(x, a n ) < (n + l/2)n. We again use the crash lemma as we did in the n = 
case to show the a„-orbit is defined for all x. The only difference is that we must show that 
9(x, a n ) > nn. That fact follows easily from lemmas 2.6b and 6.3. 




Figure 1: The projection of a-orbits on the (/, /') phase plane for several selected values of the 
shooting parameter a. 
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3 Local existence 



Proposition 3.1. There exists a two-parameter family of local solutions of Eqs.(2.14) near 
x = satisfying the initial conditions (2.17). 

Proof. The proof is standard so we just outline it. We introduce new variables w = f , z = ln(C), 
B = (1 — A)/x, and rewrite Eqs.(2.14) as the first order system 

/' - w, (3.1a) 

(x 2 w)' - ^((xf-B)w + f(l-AC 2 )), (3.1b) 

(x 2 B)' = x 2 (Aiu 2 + f 2 AC 2 + f), (3.1c) 

J = S^. (3.1d, 

We will use the sup norm throughout this discussion: \\h\\ means the sup{|/i(a;)| : < x < r}. 

Consider the space X of quadruples of functions (f,y,B,z) where ||/ — a\\ < l,||w|| < 
1, ||-B|| < M, and \\z — ln(a)|| < 1 and each of the four functions is in C°([0,r]), the space 
of continuous functions defined on the interval < x < r with the sup norm. X is a com- 
plete metric space if we take as metric the maximum of the four components. We define a map 
T : X -» X by T(f, w, B, z) = {T U T 2 , T 3 , T 4 ) where 



Ti = a + / w ds, (3.2a) 
Jo 

T 2 = ~2 I"" -j ((sf 2 - B)w + f(l - AC 2 )) ds, (3.2b) 
x Jo A 

n = 4 r s2 ( Aw2 + f 2A ° 2 + / 2 )) ds ' ( 3 - 2c ) 

x Jo 

f x 1 

T 4 = lna+ -(sf 2 -B)ds. (3.2d) 
Jo A 

One verifies easily that T does in fact take X to X and that T is a contracting map if r is 
sufficiently small, and that a fixed point of T is a solution to our equations. The proof that the 
solution depends continuously on a is also routine. □ 

4 Behavior of solutions for small a 

In this section we show that the rotation 0(x, a) of the a-orbit is arbitrarily large if a is sufficiently 
small and x is sufficiently large. 

Proposition 4.1. For any n > 0, there exists a b n such that for a < b n there is an x with 
9(x, a) > nir. 

Proof. Let / = f/a. Then, Eqs.(2.14) become 

A' = ^-—^--a 2 x(Af' 2 + AC 2 f 2 + f 2 ), (4.1b) 

C = —(A-\ + a 2 x 2 f 2 ) (4.1c) 
x A 
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with the behavior at the origin 

/(0) = 1, A(0) = 1, C(0) = a. (4.2) 

For a = (decoupling of gravity) Eqs.(4.1bc) with conditions (4.2) have constant flat-spactime 
solutions A = 1, C = a. Inserting these solutions into Eq.(4.1a) gives the Bessel equation 

(x 2 f')' + x 2 (a 2 -l)f = 0, (4.3) 

whose unique solution satisfying (4.2) is 



- sin y/a 2 - 1 x 

= r-2 — f — ' ( 4 - 4 ) 

Va —lx 

This solution has infinite rotation asi^oo. If x > then 0(x, 0) > nir so for a close to 0, 

say a < b n , we have 6(x, a) > nn because solutions of Eqs.(4.1) are continuous in a and x. This 
concludes the proof of Proposition 4.1. □ 

5 Behavior of solutions for large a 

Proposition 5.1. The a-orbits with sufficiently large a exit Q4 directly to Ql. 
We define new variables 

y = ax, v(y) = a(f(x) - a), A(y) = A(x), C(y) = C(x). (5.1) 
Then, Eqs.(2.14) become (where now the prime denotes the derivative with respect to y) 

tf) = 11(1-^(1 + J), (,2a) 
A' = l -^-y[^Ai' 2 + {l + AC-)[l + ^) \, (5.2b) 
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The initial conditions at y — are 

i(0) = 1, (7(0) = a > 1, v(0) = 0, v'(0) = 0. (5.3) 

As a — » oo, the solutions of Eqs.(5.2) tend uniformly on compact intervals to the solutions of 
the following limiting system 

yV V - y2 (1 - AC 2 ), (5.4a) 



C ) AC 

A' = l—A- v (l + AC 2 ), (5.4b) 

C = ^(A-l + y 2 ). (5.4c) 

satisfying the same initial conditions as in (5.3). The rest of this section is devoted to the analysis 
of Eqs.(5.4). Our goal is to show that v'(y) becomes positive at a point yi < y. This would 
imply that v(y) is bounded below, i.e., there is d > such that v(y) > —d for y < y\, and 
therefore v(y) > — d — 1 if a is sufficiently large. Then f(x) > a — (d + l)/a and f'(x) > for 
x = yi /a, hence, if a > y/d+l, the a-orbit exits Q4 to Ql directly without entering Q3. 

Note that the function v decouples from Eqs.(5.4bc) for the metric coefficients - this fact 
considerably simplifies the analysis. 
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Lemma 5.2. The solution of Eqs.(5.4) crashes at some y, that is, A(y) = Imiy^^ A(y) = 0. 
Moreover, 1 < y < \/3- 

Proof. Note that (yA)' < 1 — y 2 so integrating gives A < 1 — y 2 /3. Therefore, A(y) = for 

?y < v3. To show that y > 1 assume y < 1. Then (|y) = —^j- >0soif0<T<y<ywe 
have y/C(y) > r/C(r) or C(y) < C(r)/r so C is bounded. Since (AC)' = -yAC 3 , (\n(AC))' = 
—yC 2 is bounded below. Thus, by integrating one concludes that lim y ^ AC(y) > 0. But 
(AC) 2 = A(AC 2 );AC 2 < a 2 and lim^ A(y) = so \im y _ ff ( AC \y)) 2 = 0. This contradicts 
liniy^y AC(y) > so we must have y > 1. □ 

Proof of Proposition 5.1. In order to prove that i>'(y) becomes positive at some point y\ < y, 
we will show that v'(y) > 0. By Eq.(5.4a) we have v'(y) = ^ f£ z < " 1 A q C ^ dz, so we must show 

that J^ y2 ^ c2 U y>0. 

The proof of this fact is divided into two cases: (i) y 2 > 3/2, and (ii) y 2 < 3/2. Before 
considering these cases we list some useful properties of the function g = 1 — AC 2 . 

Lemma 5.3. We have: 

(a) g' = (1 - A y 2 g)C 2 /y; 

(b) if g(yo) > 0, then g(y) > for all y > y ; 

(c) g'>0ifg<l/3; 

Proof. Part (a) is a calculation. For (b) note that g'(g = 0) > so g cannot cross zero from above. 
For (c) we have (yA) 1 < 1—y 2 so integrating gives 1 — A > y 2 /3 and hence, g' > y(l/3—g)C 2 . □ 

We return now to the proof that ^j§dy > 0. We first consider the case (i) y 2 > 3/2. 

A calculation shows that y 2 g = (2y 3 /3 + yA- y)', hence y 2 gdy = 2y 3 /3-y > if y 2 > 3/2. 
Since g(0) = 1 — a 2 < 0, this implies that g(a) = for some a < y and therefore g(y) > for 
y > a. Note that AC is monotone decreasing because (AC)' = —yAC 3 < 0. Thus 

V 2 9(y) for < y <y (5 . 5) 



A(y)C(y) ~ A(a)C(a 

and therefore 



□ 

Now we consider the case (ii) y 2 < 3/2. 

Lemma 5.4. Define the function p = 1 + y 2 g — y 2 ■ If y 2 < 3/2, then p(y) > 0. 

Proof. Note that p(0) = 1. Let j/i be the first zero of p, that is, p(yi) = and p'(yi) < 0. If 
g(yi) > 1/3 then p = y 2 g + 1 — y 2 > y 2 /3 + 1 — y 2 = 1 — 2y 2 /3. Thus p can have a zero for 
y\ < 3/2 only if g(yi) < 1/3. Then, from Lemma (5.3), g'(y) > for all y < y\. Define a 
function k(y) = 2 — 3A — y 2 . A calculation gives y 3 g' = (y(k + p))' so by integrating we get 
k(yi) > 0. On the other hand we have p' = yC 2 (k —p), so k(y\) < 0; contradiction. □ 

To show that ^j§dy > 0, we rewrite it as 

fy y 2 g f y P -l + y 2 , f y p , fl-y 2 , 
Jo AC"" = Jo —^ dV = Jo AC'' ~ Jo (5J) 
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The first term on the right hand side of (5.7) is positive because p is positive. To compute the 
second term, note that 



hence L = lim y ^y(y/C) exists and is finite since y > 1 by Lemma 5.2. If L > then 
linij,^ C(?y) = y/£ < oo so Cis bounded. Since lim y ^^ A(y) = we conclude that lining AC{y) 
0. But (In AC)' = yC 2 is bounded so In AC is bounded below and hence lim AC ^ 0. This con- 
tradiction shows that L = 0. Thus, the second term on the right hand side of (5.7) is zero. This 
concludes the proof of Proposition 5.1. 



6 Asymptotics of globally regular solutions 

In this section we derive the leading asymptotic behavior of globally regular solutions. We use 
lim to denote lim^^oo. 

Proposition 6.1. An a-orbit which exists for all x and has bounded rotation is asymptotically 
flat. The leading asymptotic behavior for x — > oo is 

A (x)~l-?^, C{x)^C 00 e^, f(x) ~ Ue-**, (6.1) 

where < < oo, < C x < I, and b = y/1 - C£_, . 

To prove this proposition we need several partial results. 

Lemma 6.2. An a-orbit which exists for all x and has bounded rotation is ultimately in the 
second (Q2) or fourth (Q4) quadrant. 

Proof. If 9{x) is bounded above then there is an integer n > such that 9(x) < (n + l/2)n 
for all x but 6{x\) > (n — l/2)it for some x\ and hence, by Lemma 2.6a for all x > x\, 
(n — l/2)ir < 9(x) < (n + l/2)n. Wc next show that there is an X2 such that for all x > X2, 
nir < 6(x) < (n+l/2)n (that is, the orbit is ultimately in Q2 or Q4). Note that, by Lemma 2.6c 
the orbit must satisfy either nir < 6(x) < (n + l/2)n or (n — l/2)ir < 6(x) < mr, that is the 
orbit must lie in Q3 or Q2 if n is odd and in Ql or Q4 if n is even. We must rule out the 
possibility that the orbit is in Ql or Q3. Assume that the orbit lies in Ql or Q3 for all x > x\. 
Then f(x)f'{x) > for all x > xi, so f 2 (x) > f 2 (xi) for all x > x\. From Eq.(2.15b) we have 
(xA)' = 1 - x 2 Af' 2 - x 2 f 2 AC 2 - x 2 f 2 so (xA)' < 1 - x 2 f 2 < 1 - x 2 f 2 ( Xl ) and hence A goes 
to zero in finite x. This contradiction concludes the proof. □ 



Lemma 6.3. Under the assumptions of Proposition 6.1 the function g = 1 — AC 2 is eventually 
positive. 

Proof. Suppose that g(x) < for all x. We claim that this implies lim A =1. To see this, 
suppose that liminf A = 1 — 4e for some e > 0. Let — /3 — lim g < which exists because g' > 0. 
Note that g(x) < —[3 for all x. Choose an X\ such that g{x\) > —(3 — e. If A(x 2 ) < 1 — 3e 
for some x 2 > x 1 , then by (2.18) g'(x) > C 2 (l - A)/x > (1 - A)/x = x(l - A)/x 2 > x 2 (l - 
A(x 2 ))/x 2 > 3ex 2 /x 2 for x > x 2 , where the last but one inequality follows from the fact that 
x(l — A(x)) is monotone increasing. Integrating this inequality from x 2 to 2x 2 say, we get 
g(2x 2 ) > g(x 2 ) + 3e/2 > — [3 — e + 3e/2 > —(3; contradiction. Thus, liminf A = 1 and hence 
limA = 1. Since limg = lim(l— AC 2 ) exists, lim C also exists and is finite. Next, fromLemma6.2 
we know that the a-orbit is ultimately in Q2 or in Q4. For concreteness we consider the case of 
Q4 (the proof of the Q2-case is identical), that is f{x) > and f'(x) < for sufficiently large 
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x. Then, from (2.14a), lim(x 2 /'/C) exists so lim(x 2 /') = — t < exists as well (where r might 
be infinite; the point is that r 7^ 0). Now, by L'Hopital's rule, lima;/ = — lim(x 2 /') = r. But 
(2.14c) says (InC)' > t 2 /4x which implies limC — 00, a contradiction. □ 

Proof of Proposition 6.2. From the previous lemma we know that there exists an x\ such that 
g{x) > for x > x\. Let u = ACf/g for x > x\. A calculation shows that u' = —AC(fC 2 (l — 
A)/x - f'g + xff 2 )/g 2 so u' < if g > 0. Multiplying Eq.(2.14a) by u we obtain 

(x 2 Aff'/ 9 y = x 2 f 2 +x 2 f'u'/C. (6.2) 

The right hand side is positive for x > x\ so x 2 Aff'/ g is negative and increasing, hence it has a 
finite non-positive limit. This implies that x 2 f 2 is integrablc. Similarly, multiplying Eq.(2.14a) 
by / we obtain 

x 2 ff'/C = (x 2 f 2 g + Ax 2 f' 2 )/(AC). (6.3) 

The right hand side is positive for x > x x so x 2 ff'/C is negative and increasing, hence it has 
a finite non-positive limit. This implies that Ax 2 f' 2 is integrable (recall that AC is monotone 
decreasing). The integrability of x 2 f 2 and Ax 2 f implies via Eq.(2.15) that limM = Moo < 00 
exists. This concludes the proof that A(x) <~ 1 — 2Moo/x. 

Having limyl = 1 we can strengthen Lemma 6.3 by showing that lim g = <?oo > exists. To 
see this choose an x\ such that g{x\) > 0. Then AC 2 (xi) < 1, hence AC(xi) < 1. Since AC 
is monotone decreasing, we have AC(x) < AC(xi) for x > xi and thus lim^4C < 1. Hence, 
lim^4C 2 = (lim AC) 2 / 'lim A < 1. Since g =1 - AC 2 , limg exists and limg > 0. 

Now we have all we need to derive the asymptotics of /. Let r = f'/f. Then r' = f" / f — r 2 = 
— r(l + A — x 2 f 2 )/(xA) + g/A = g^ — r 2 + e(x), where lime = 0. Let a(x 2 ) = max(|e(a;)|) for 
x > X2 and assume that X2 is sufficiently large so that g^ > cr(x2). If r(x2) > — \jg0o~0\x2), 
then clearly r becomes eventually positive which contradicts that the orbit is eventually in Q2 
or Q4. If r(x 2 ) < — \fgZ> + 0^2), then limr = —00 - this is impossible because then by 
L'Hopital's rule limr = lim/"//' = limg/r = 0. Therefore r(x2) must be sandwiched in the 
interval —\fg~Z) + o~{x2) < ^(^2) < —\/goo — o-{x2). Since X2 is arbitrarily large and limer = 0, 
we conclude that limr = —^goo- The asymptotics of / given in (6.1) follows immediately from 
this. 

Finally, inserting the derived leading asymptotic behavior of A and / into Eq.(2.14c), we 
obtain C'/C <~ —IMoejx, from which the asymptotics of C follows trivially. □ 

7 Solutions that crash 

Proposition 7.1. // the b-orbit crashes at some x then g(x) > for x near x. 

Proof. Suppose that g{x) < for all x < x, so AC 2 (x) > 1 for all x < x. We have from (2.18) 
that g' > AC 2 xf > xf . Integrating this inequality from some X\ > to some x 2 < x, we 
obtain 

rx 2 rx-2 

xi / fdx < / xf' 2 dx < g(x 2 ) - g(xi) < a 2 - 1, (7.1) 

which implies (by the Cauchy-Schwartz inequality) that / is bounded. 

Next, A(x) — 0, AC 2 > 1, implies that \im x ^ s - C = 00; moreover, by (2.14c) (InC)' < 
xf 2 /A, hence xf 2 /A is not integrable near x. Since / is bounded, this shows that 1/A is not 
integrable near x. But from (2.18), g' > C 2 (l - A)/x = AC 2 (1 - A)/(xA) > l/(2xA), so g' is 
not integrable near x, which contradicts the fact that g is a bounded function. □ 
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The importance of Proposition 7.1 derives from Lemma 2.6b which says that if g > then 
rotation stops. The main result of this section is the crash theorem which states that if an orbit 
has bounded rotation and crashes, then nearby orbits also have similarly bounded rotation. The 
precise statement is given in Proposition 7.2. Since we consider more than one orbit in this 
section, we use the notation A(x, a) to denote the value of A at x for the a-orbit, etc. 

Proposition 7.2 (Crash Theorem). If the b-orbit crashes at x — x and 

(a) if (k — 1/2)tv < 6(x, b) < kn, k > 1, for x near x, then nearby orbits have rotation < kn for 

x > x; 

(b) if kn < 9(x,b) < (k + l/2)n, then nearby orbits have rotation < (k + l/2)n. 

Proof. Part (a): Suppose the 6-orbit crashes in Q3 or Ql. By Proposition 7.1, g(x\,b) > for 
some X\ < x with (k— l/2)7r < 9(x\, b) < kn; hence, for a sufficiently near b we have g(x\, a) > 
with (k — l/2)n < 6{x,a) < kn. By Lemma 2.6b, 0(x,a) < kir for all x > X\. Part (b): This 
case is much more difficult and will require several auxiliary results. It follows from part (a) 
that nearby orbits have rotation < (k + l)ir; we must prove a much more difficult result, namely 
that nearby orbits have rotation < (k + l/2)n. □ 

Remark. It is clear from numerical observations that no a-orbit crashes in Q2 or Q4; however, 
that appears to be quite difficult to prove. Moreover, one can easily construct orbit segments 
that start, for example, at x = 1 with / = 5, /' = 0, A = 0.2, C = 3, say, that crash in Q4. 
Such orbit segments have \im x ^ x - f'(x) = — oo. Nevertheless, the next lemma shows that Af' 2 
remains bounded at crash. 

Lemma 7.3. If an a-orbit is defined for x < X2,ff'(x) < for X\ < x < x 2 , f 2 (xi) < B, 
and f'(x\) — then Af' 2 (x) < max(B, a 2 /3). In particular, if an orbit crashes in Q2 or Q4, 
lirn^s- A(x)f'(x) = 0. 

Proof. We set q = Af' 2 and then compute that 

xq' = -(3 + x 2 /' 2 + x 2 C 2 .f 2 )q f + 2xff> + x 2 f f 2AC 2 xff. (7.2) 

Note that q > and all terms on the right side of (7.3) are negative except for the last two. 
If q > B, we combine the term —qx 2 f' 2 with x 2 f 2 f' 2 ; clearly, x 2 f 2 f' 2 — qx 2 f' 2 = (f 2 — 
q)x 2 f' 2 < 0. Next, we combine the term -qx 2 f' 2 C 2 with —2xff'AC 2 to get -AC 2 (y 2 - 2y) 
where y = —xff; the maximum value of this expression occurs when y = 1 and that value 
is AC 2 < a 2 by Lemma 2.4. Hence, if q > a 2 /3, then -q(x 2 f' 2 C 2 ) - 3q - 2xff 2 AC 2 < 0. 
Thus, q > max(£>, a 2 /3) implies that q' < 0; consequently, Af' 2 (x) < max(_B, a 2 /3). Since 
AAf' 2 = (Af) 2 , and Af' 2 is bounded and limx ^ x-A(x) = 0, lim x ^ 2 - (A(x)f'(x)) 2 = 0, 
hence \im x ^ s - A(x)f'(x) =0. □ 

We can now discuss the strategy of the proof of part (b) of Proposition 7.2. We want to show 
that if an orbit is sufficiently close to an orbit that crashes in Q4 then it must either crash or 
exit Q4 to Ql (the case in which the orbit crashes in Q2 is completely symmetric). To that 
end, let v(x) = A(x)f'(x). We will prove that v(x, a) goes to if a is sufficiently close to b and 
f(x, a) > 0. This means either /' = and hence the orbit is exiting Q4 to Ql, or A = 0, that is, 
the orbit is crashing in Q4. Note that v'(x) = -(2Af'-xf+xAC 2 f + x 2 Af' 2 +x 2 f 2 f'AC 2 )/x = 
—v(2 + x 2 f' 2 + x 2 f 2 C 2 )/x + fg > fg. We know that v(x, b) goes to at crash so nearby orbits 
will also have v small for x near x. We will show that / and g are both uniformly bounded away 
from in an interval about x. That is, the size of the interval and the bounds work for all a near 
b. That is enough to force v positive. The most technical part of the proof involves showing that 
nearby orbits stay in Q4 long enough to have v go positive. Since /' goes to — oo at crash, nearby 
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orbits have /' large also. Now, (2.14a) can be written as xAf" + (1 + A — x 2 f 2 )f — xgf = 0; 
moreover, to get to Q3 orbits must pass through xf(x) < 1 which means that the coefficient of 
/', (1 + A — x 2 f 2 ), is positive. That is enough to bound /'. 

The details of the proof, especially Lemma 7.5, are tedious. We will restrict ourselves to an 
interval 0.99 x < x < 1.01 x and replace x by x (whenever justified) in making estimates. 

We show next that if the o-orbit crashes at x — x with rotation kn < 6(x, b) < (k + l/2)w 
then \xf(x)\ > 1. 

Lemma 7.4. // the b-orbit crashes at x — x with 9{x,b) < (k + l/2)n for all x < x and 

9(x,b) > kn for x near x, then \xf{x)\ > 1, in particular f(x) ^ 0. 

Proof. The assumption on Q(x, b) tells us that the orbit lies in Q2 or Q4 for x near x. For 
simplicity of exposition we only discuss the case of Q4, i.e., f(x) > 0,f'(x) < 0. In particular, 
/ is a monotone function and hence has a limit at x. Thus, h{x) = xf{x) is continuous; in 
particular, if we suppose that xf{x) < 1, then h(x) < 1 for x near x. Since A{x) = 0, we get 
from (2.14c) that xAC = C(A — 1 + x 2 f 2 ) < for x near x. We conclude that C is bounded 
above, hence lim^^g- AC 2 = and lim^^g- g = 1. Since g > 0, the right hand side of Eq. (2. 14a) 
is positive and hence x 2 f'/C is bounded and since C is bounded we conclude that /' is bounded; 
thus lim^- Af' 2 = 0. Then, from (2.14b), xA' = 1 - A - x 2 f 2 - x 2 (Af 2 + AC 2 f 2 ), we see 
that A' > near x so there is no crash. This is a contradiction so we conclude that xf(x) > 1 
and hence f(x) > 0. □ 

Lemma 7.5. There is a 7 > such that h{x,a) — xf(x,a) > 1/4 for all a sufficiently near b 
and x < x < x + 7. 

Proof. If the &-orbit crashes at x = x with rotation 8(x, b) > kn, then there is a y such that 
9(y,b) — kn. Let B = (f(y,b) + l) 2 . By Proposition 7.3, if a is sufficiently close to 6, Af 
(along the a-orbit) is bounded in Q4 by D = max(a 2 /3, £>); D is a uniform bound on Af' 2 
in Q4 for all a sufficiently near b. Next, choose x\ such that 0.99 x < x\ < x and such 
that A(xi,b) < 0.01, g{x\, b) — 2r > 0, and h(xi,b) > 0.9; this is possible by Lemma 7.4 
and Proposition 7.1. Then, for a sufficiently near b we have A(xi,a) < 0.02, g(x\, a) > t > 
0, f(xi,a) < f{x 1 ,b) + 0.01/5 and h{x x ) > 3/4. We shall find a 7 G (0,0.015) that works for 
all a, that is, it satisfies h(x,a) > 1/4 for all a sufficiently near b and x < x < x + 7. So let 
a satisfy: i) Af' 2 (along the a-orbit) is bounded by D, ii) A(xi,a) < 0.02, iii) h(x\,a) > 3/4, 
and iv) g(xi, a) > t > 0. If /i(x, a) > 1/4 for all x < 1.01 x and all a near b we are done - let 
7 = 0.01 x. Otherwise, we define x-i = x-iia), etc. by /i(x 2 ) = 3/4, h(x 3 ) = l/2,/i(x 4 ) = 1/4, 
where x 2 ,X3, and X4 are the largest values of x < 1.01 x with that property. For x > x 2 we 
have from (2.14a) xAf" = xgf - (1 + A - h 2 )f > -(1 + A — h 2 )f > - f/4 since h < 3/4 so 
/" > -f' 3 /(4xAf 2 ) > -/ /3 /(4 * 1.01 xD) or f"/f' 2 > -/'/(4.04 x D). We now integrate the 
above from x-i to x > x-$ to get 



/'(a) - f'{x) + f'(x 2 ) ~ J X2 f> 2 X ~ J X2 AMxD aX 4.04 x D ' 



Now, f(x) > f(x 3 ), so -f'(x) < f(x2 f_ D f(x3) w h(x %h (x3) = 20x 2 C. Using the uniform bound 
on /' in the interval x 3 < x < X4, we have X4 — x 3 = {f{xi) — /(£3))//'(£) for some ^ e [x%, X4}. 
But (/(x 4 ) - /(z 3 ))//'(0 > (h(x 4 ) - h(x 3 ))/(xf(t)) > 1/80 x 3 D and hence we may take 
7 = l/80x 3 L>. □ 

Lemma 7.6. In the interval X\ < x < x + 7, g(x,a) > min(r, 0.9/h 2 (x, bj). 
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Proof. From (2.18) we have xg' = C* 2 (l -A + x 2 Af 2 - x 2 gf 2 ) > C 2 (l- A- x 2 gf 2 ). Moreover, 
since A{x\,a) < 0.02 and xA' < l,A(x,a) = A(x 1 ,a) + A'(z)(x-x 1 ) < 0.02 + l/z(0.02x) < 0.04, 
so if g < 0.96/h 2 (x) then g' > 0. Since f(xi,a) < f(xi,b) +0.01/x,h(x,a) < 1.01xf(xi,a) < 
1.01(xf{xi,b) + 0.01) < 1.02 xf(xi,b), we have g' > if g(xi,a) < 0.9/h 2 {x,b). Thus, if 
t < g(xi,a) < 0.9/h(x,b), g' > 0, and g(x,a) > r in the interval x\ < x < x + 7; if g(xi,a) > 
0.9/h 2 (x,b), then g(x,a) > 0.9/h 2 (x,b) for all x in the interval x\ < x < 2 + 7 because g cannot 
cross that value from above. □ 

Note that the above lower bound on g is uniform - it applies to all a satisfying the conditions 
i) Af' 2 (along the a-orbit) is bounded by D, ii) A(xi,a) < 0.02, iii) h(x\,a) > 3/4, and iv) 
g(xi,a) > t > 0. 

Lemma 7.7. For all a sufficiently near b, v(x, a) goes to for some x < x + 7. 

Proof. To show that v(x, a) goes to 0, we note that h(x,a) > 1/4 for all a near b and x < x < x+j 
by Lemma 7.5. Hence, f(x,a) = h(x,a)/x > 1/Ax. By Lemma 7.6, g(x,a) > min(r, l/h(x), 
hence v' > l/4Smin(r, l/h(x) = n > for x < x < x + 7. Thus, v(x + 7) — v(x) = v'dx > 
/? +7 Tjdx > 777. Let X\ be chosen so that v(xi,a) > —^7/2. Then, if a is sufficiently close to b, 
\v(xi, a) — v(xi, b)\ > 777/2 so v(xi, a) > —777. For such a we then have v(x + -f, a) > v(x, a) + 777 
and 7j(S, a) > v(x\,a) > —777 because v' > fg > 0; thus, v(x + 7, a) > 0. □ 

We now complete the proof of Proposition 7.2. 

Proof of Proposition 7.2 b). Suppose that the &-orbit crashes at x — x with 6(x, b) < (k + 1/2)tt 
for all x < x and 6(x, b) > kir for x near x. For a near b there is an x < x + 7 with u(x, a) = 
by Lemma 7.7. Since a; < x + 7, h(x) > 1/4, i.e., f(x,a) > 0, so the a-orbit crashes, A(x,a) = 0, 
or exits Q4 to Ql (or Q2 to Q3), f'(x,a) = 0, never to return. In either case, the a-orbit has 
rotation 9(x, a) < (k + l/2)n. □ 



8 Proof of the main theorem 

Proof of Theorem 2.1. Let X n = {a > | 9(x, a) < (n + l/2)n for all x for which the a-orbit 
is defined}. Note that X n _\ C X n and X n 7^ by Proposition 5.1 and hence, X n ^ 0. Also 
note that b n+ \ > is a lower bound for X„ by Proposition 4.1; hence, X„ has a greatest lower 
bound a„ = mf(X n ) > b n+ i > 0. We will show that the a„-orbit is a globally regular solution 
and 7i7r < 9(x, a n ) < (n + l/2)7r for large x. 

We first show that a n G X n , i.e., a n is the smallest element in X n . If 0(x, a n ) > (n + 1/2)tt 
for some x then 0(x, a) > (n+ l/2)n for all a near «„ so « ^ X„ for these a's and this contradicts 
the fact that a n is the greatest lower bound of X n . Thus, a n G X„. In particular, the a„-orbit 
has bounded rotation. 

Next we show that the a„-orbit does not crash. Recall from Proposition 7.1 that if the a n - 
orbit crashes at x = x then g(x,a n ) > for x near x. If the a„-orbit crashes in Ql or Q3, that 
is, if 9{x, a n ) < mr for x near x then 9(x, a) < nir and g(x, a) > for all a near a n which implies 
by Lemma 2.6b that the a-orbit must have 9(x, a) < mr for all x. Thus, a G X n for all a near 
a n and this contradicts the fact that a n is the greatest lower bound of X n . 

Similarly, if the a„-orbit crashes in Q2 or Q4, that is, at some x with (n+ l/2)7r > 9{x, a n ) > 
n7r, then by the crash lemma (n + 1/2)tt > 8(x, a) for all x in the domain of definition of the 
a-orbit for all a near a n and this contradicts the fact that a n is the greatest lower bound of X n . 

Thus, the a„-orbit is defined for all x and hence is a globally regular solution by Proposi- 
tions 6.1. Also, by Proposition 6.2, the a„-orbit is in Q2 or Q4 for large x. It remains to prove 
that 9(x,a n ) > nn for large x. Suppose that 9(x,a n ) < mr for large x. By Lemma 6.3 we have 
that g(x, a n ) > for large x and hence, g(x, a) > for all a near a n . Then, by Lemma 2.6b the 
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a-orbit must have 6{x, a) < nir for all x and thus a G X n , and this contradicts the fact that a n 
is the greatest lower bound for X n . This completes the proof of Theorem 2.1. □ 
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